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Abstract 

We consider a transmission wave equation in two embedded domains in R^, where 
the speed is ai > in the inner domain and a2 > in the outer domain. We prove a 
global Carleman inequality for this problem under the hypothesis that the inner domain 
is strictly convex and oi > 02. As a consequence of this inequality, uniqueness and Lip- 
schitz stability are obtained for the inverse problem of retrieving a stationary potential 
for the wave equation with Dirichlet data and discontinuous principal coefficient from a 
single time-dependent Neumann boundary measurement. 
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1 Introduction and main results 
1.1 Presentation of the problem 

The inverse problem of recovering coefficients from a wave equation with discontinuous co- 
efficients from boundary measurements arises naturally in geophysics and more precisely, in 
seismic prospection of Earth inner layers [20] . 

Here we are interested in the case where only one particular measurement is available. 
This could be important, for instance, in seismic prospection, where data of a single wave 
that propagates through the Earth is considered. 

Consider two embedded domains, where the speed coefficients are ai > in the inner 
domain and 02 > in the outer domain. Stability of the inverse problem we study here is 
obtained by deriving a global Carleman estimate for the wave equation with discontinuous 
coefficients. We prove this Carleman inequality in the case the inner domain is strictly convex 
and the speed is monotonically increasing from the outer to the inner layers, i.e. oi > 02. 
This last situation is, incidentally, the general case into the Earth. 

Figure 11.11 illustrates the role of these hypothesis and gives some intuition with the help 
of Snell's law. In the case oi > a2 (see Figure II. ![ left) the incident rays coming from the 
inner domain toward the outer domain become closer to the normal at the interface since 
sin(0i) > sin(6'2), where 9i, i — 1,2 are the corresponding incident angles. Therefore, all 
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the rays coming from the inner baU with any incident angle 9i in (— 7r/2,7r/2) succeed in 
crossing the interface. In the opposite case, when ai < 02 (see Figure ll-H center) we have 
sin(6'i) < sin(6'2) and there is a critical incident angle 9* < 7r/2 such that the rays with 
incident angles 9i out of the range {—9*, 9*) remain supported near the interface and do not 
reach the outer domain, so this information does not arrive at the exterior boundary. Finally, 
strict convexity of the inner domain avoids trapped rays (see Figure [TTTl right). 




Figure 1: Two domains with speed coefficients ai (inner) and a2 (outer). In the first figure (left), if 
ai > a2 by Snail's law all the inner rays reach the exterior boundary independently of their incident 
angles. Conversely, in the second figure (center) if ai < a2 some rays with large incident angles 
remain trapped near the inner interface. The last figure (right) shows a trapped ray into a captive 
domain. 

Global Carleman estimates and the method of Bukhgeim-Klibanov [9] , [7] are especially 
useful for solving the one measurement inverse problems. It is possible to obtain local Lips- 
chitz stability around the single known solution, provided that this solution is regular enough 
and contains enough information J^7\ (see also and [Hj). Many other related inverse re- 
sults for hyperbolic equations use the same strategy. A complete list is too long to be given 
here. To cite some of them see [36] and [41] where Dirichlet boundary data and Neumann 
measurements are considered and |23j . |24j where Neumann boundary data and Dirichlet 
measurements are studied. These references are all based upon the use of local or global 
Carleman estimates. Related to this, there are also general pointwise Carleman estimates 
that are also useful in similar inverse problems [THl dZl [55] ■ 

Recently, global Carleman estimates and applications to one-measurement inverse prob- 
lems were obtained in the case of variable but still regular coefficients, see [22] for the isotropic 
case, and [3D] and [3] for the anisotropic case. It is interesting to note that these authors 
require a bound on the gradient of the coefficients, so that the idea of approximating discon- 
tinuous coefficients by smooth ones is not useful. 

There are a number of important works [371 1191 1381 1391 [5] concerning the same inverse 
problem in the case that several boundary measurements are available. In these cases, it is 
possible to retrieve speed coefficients and even discontinuity interfaces without any restrictive 
hypothesis of strict convexity or speed monotonicity. For instance, one can retrieve the 
interface by observing the traveltime reflection of several waves. Indeed, it is well known 
that the interface can be recovered as the envelope of certain curves as shown in Figure 11.11 
(see also [20] and the references therein). This method works independently of the sign of 
oi — 02 and this explains in part why there are no geometrical or speed monotonic hypotheses 
for these kind of inverse results. 

Let us now give some insight into the relationship between this work and exact control- 
lability or energy decay for the wave equation with discontinuous coefficients. 

First of all, the global Carleman estimate we obtain immediately implies a particular 
case of a well known result of exact controllability for the transmission wave equation [31] . 
Roughly speaking, the result of |^ states that we can control internal waves from the 
exterior boundary in a layered speed media if the speed is monotonically increasing from the 
outer to the inner layers and the inner domain is star shaped, a weaker assumption than 
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Figure 2: Recovering the interface as the envelope of circumferences. Each circumference, centered 
at some point x on the exterior boundary, represents the possible locations of the nearest point of 
the inner interface where the reflection took place for a given traveltime measured at x. 



strict convexity. Moreover, if the speed monotonicity is inverted, there are non controllable 
solutions with concentrated energy near the interface ^3]; [12] . 

Secondly, there exist several results about the growth of the resolvent for the spectral 
stationary transmission problem, from where it is possible to derive the speed of local energy 
decay for the evolution wave equation with transmission conditions [lOj . In the case oi > 02 
and if the inner domain is strictly convex, it has been shown using micro-local analysis 
[34l [35] that the speed of the energy decay is exponential if the dimension of space is odd 
and polynomial otherwise. In the general case, including the cases when oi < 02 or the 
inner domain is not strictly convex, it has been proved using micro-local analysis and global 
Carleman estimates for the spectral problem [3] that the energy decays as the inverse of the 
logarithm of the time. 

Notice that we shall only consider here the case of a discontinuous coefficient which 
is constant on each subdomain (i.e. ai and 02 constants). We will indeed concentrate our 
discussion on the main difficulty, namely the discontinuity at the interface. However, we could 
also consider variable coefficients ai{x) and 02(0:) such that their traces at the interface are 
constant, under additional assumptions of boundedness of Voj similar to those appearing in 
[22] (see Remark |4|). 

Finnaly, we note that a global Carleman estimate [15] has also been obtained for the heat 
equation with discontinuous coefficients. That work was initially motivated by the study 
of the exact null controllability of the semilinear wave equation, but the estimate has been 
recently used to prove local Lipschitz stability for a one measurement inverse problem for 
the heat equation with discontinuous coefficients [5] , [5] . 

Having introduced the problem, let us now present our main results. 

1.2 Inverse problem 

Let Vt and ili be two open subsets of with smooth boundaries F and Fi. Suppose that 
fii is simply connected, fii C and set 5^2 = \ ili. Thus, we have dVl2 — F U Fi. We also 
set: 




with flj > for j = 1, 2. We consider the following wave equation: 




(x, t) e r2 X (0, T) 
(x, t) e F X (0, T) 
uq a; e 17 



u 



(1) 



u(0) 
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We know that [32l [13] for each p g L°°(17), uq G HQ{fl) and ui e i^(il), there exists a 
unique weak solution u{p) of equation ^ such that 

u{p) e C([0,T];i/i(r!)), ut e C{[0,T];L\n)). 

We shall prove the well-posedness of the inverse problem consisting of retrieving the potential 
p involved in equation ([1]), by knowing the flux (the normal derivative) of the solution u{p) 
of ([1]) on the boundary. We will prove uniqueness and stability of the non linear inverse 
problem characterized by the non linear application 



, du 



(2) 

rx(o,T) 



More precisely, we will answer the following questions. 
Uniqueness : 

Does the equality = on T x (0, T) imply q = p on Q 7 



Stability 



Is it possible to estimate (q — by 



(q) _ du(p) 



in suitable norms ? 

rx(o,T) 



The idea is to reduce the nonlinear inverse problem to some perturbed inverse problem 
which will be solved with the help of a global Carleman estimate. More precisely, we will 
give a local answer about the determination of p, working first on the perturbed version 
of the problem, as shown is Section 3. Assuming that p G L°° is a given function, we are 
concerned with the stability around p. That is to say, p and u{p) are known while q and u{q) 
are unknown. 

We are able to prove the following result, which states the stability of the inverse problem. 

Theorem 1 Assume is bounded, fli CC fl is a strictly convex domain with boundary Ti 
of class and ai > a2 > 0. There exists Tq > such that, given T > Tq, if p G L°°{fl), 
Uq £ iJQ(il), ui G L'^{fl) and r > satisfy 

• |uo(ic)| > r > a. e. in fl, and 

• uip) e H\0,T;L°^{n)) 

then, given a bounded setlA C L°°{^T}, there exists a constant 

C = C(ai,a2,f^i,f^2,T, |b||L°°(n), \\u{p)\\ m (^l^) ,U , r) > such that: 



du{p) du{q) 



02^- 02 



diy dv 



for all q £lA, where u{p) and u{q) are the solutions of (QP with potential p and q, respectively. 

Remark 1 In section\^ is given an estimate for Tq in function of ai, 02, and ^2- See 
Theorem O 

Let us remark that as a direct consequence of the local stability of Theorem [T] we have 
the following global uniqueness for our inverse problem: 

Corollary 1 If u{p) and u{q) are two solutions of {IJ) for potentials p and q in i°°(r2) with 
u{p), Uq, Ui, ai, 02, fli, fl2 and T satisfying the hypothesis of Theorem]^ and such that 
9^ = ^ on Fx (0, T) then p = q. 

The proof of this result consists of two parts: a global Carleman estimate and the reso- 
lution of the inverse problem and lipschitz stability following the methods introduced in [9] 
and [27] which we have already mentioned in the introduction. 
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1.3 Carleman estimate 

We introduce here our main result concerning a global Carleman estimate for the solutions of 
problem ^ extended to the time interval (-T, T). We set g = 17 x (-T, T), S = T x (-T, T), 
Si = Fi X [~T,T), Qj = Qj X [—T,T), Uj =^ wig and i^j the outward unit normal to Qj, 
for j = 1,2. 

We will work with an equivalent formulation of Notice that for each / e L^{Q), u 
solves the equation 

utt — div{aS/u) + pu — f in Q (3) 
if and only if, for each j E {1, 2}, uj solves (see |31) ) 

Uj,tt - ajAuj +puj = ftq. in Qj (4) 

together with the transmission conditions 

ui = U2 on El 

a^*ii +021^^0 on El. 

In order to construct a convenient weight function, take xq € fii and for each x G r2\ {xo} 
define i{xo, x) = {xq + X{x — xq) : A > 0}. Since fli is convex there is exactly one point y{x) 
such that 

y{x) eTineixo,x). (6) 
We define the function p : \ {a;o} — » by: 

p{x) = \xo-y{x)\. (7) 

Let e > be such that C ili (and small enough in a sense that we will precise later) and 
let < El < £2 < Then we consider a cut-off function 77 e C°°(M) such that 

0<?7<1, 7] = 0m B,,{xo), r,^lmn\B,,{xQ). (8) 

For each j £ {1, 2} wc take k such that {j, k} — {1, 2} and we define the following functions 
in the whole domain x M 

(l3j{x,t)^7]{x)-^\x-xa\^-(3t^ + Mj {x,t)enxR, (9) 
p{xy 

where (3, Mi and M2 are positive numbers that will be chosen later. Then, the weight 
function we will use in this work is 

r (l)i{x,t) (x,t) e f^i X R 

Hx,t) = { (10) 

[ (t)2{x,t) {X,t) e 1^2 X R. 

Notice that (see (c) and (d) in Proposition [T] below) (f>i and (f>2 satisfy ^ if and only if 

Ml - M2 = ai - 02. (11) 

We denote 

L = d^-aA and £;(z) = |ztp - a|Vzp. 
As usual, we do the change of variables 

= e^"^ , A > 0, w = e^^u , s > 0, P(u;) = e"^L(e-'"^w) (12) 

and after algebraical computations, we split P{w) into three terms as follows: 

P{w) = Pi(w) + P2(w) + Riw), 
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where for some fixed real number 7 e (0, 1) 

Pi{w) = wtt - aAw + s'^X'^ip'^E{(l))w, 

P2{w) = {'^-l)s\ipL{4>)w- s\^ipE{4>)w-2sXif{4>tWt-aS/(j)-yw), 
R(w) = —jsXipL{(j>)w. 

We will write P'^, Pf, P^ , etc. if we want to make the dependence on (p explicit. Also, given 
U C M^, we define the norm in H^{U x {-T,T)), given by 

(•T 



T 

-T JU 



i\gt{' + \V9\'')^ + s'X' 



T JU 



and 



(13) 
(14) 



j:1 = {(x, t) e S : V^(a;, t) ■ iy{x) > 0}. 
Finally, we define the space 

X = {ueL\-T,T;L\n)) \ Luj e L\-T,T; L^Qj)), j ^ 1,2; u\s = 0, 

u{±T) = Ut{±T) = 0, and u satisfies ©}. 

The main global Carleman estimate is the following 

Theorem 2 Assume fli is a strictly convex domain of class , and ai > 02 > 0. Let 
Xk & ^1, k — 1,2 and let (f)^ , ip^ , he the corresponding functions defined for Xk as we did 
before for Xq m l\10^ and (jjj^p . Let v be the unit outward normal to il. Then there 

exists C > 0, Sg > and Aq > such that 
2 



2 



fc=i ^ 



P^ {w^) 



k=l 



LHQ) 



+ sX 



dw 



02- 



dv 



(15) 



for all u X , A > Aq and s > sq. 

Notice that in the right-hand side of we have the term 



pi' [w'^) 



LHQ) 



Since we consider an equation given by the operator Lp — du — aA + p, it is important to 



e^'^'^ \Lpu\ is also true for 



note that the same estimate with right-hand side equal to 
all potentials p such that |p|Loo(n) ^ fT^^ with m already fixed. Indeed, 

iLwp < 2\Lpu\'^ + 2m'^\u\^ 
and taking s large enough, the left hand side of the Carleman estimate of Theorem [3] can 
absorb the term 2Cm^ / / e^^'^jup. That is, we have the following result. 



Corollary 2 Under the hypothesis and notations of Theorem\^ given m G R, there exists 
C > (depending on m), sq > and Aq > such that for all p G L°°{Q) with WpllL-^^n) ^ 
we have 

2 



fc=i 



2 



< 



fc=i 



P2^ (u;'=) 
e^'"^''\Lpu\'^ + sX 



2 



dw 



02- 



dv 



(16) 



for all u G X, A > Aq and s > Sq. 
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The paper is organized as follows. In Section [2] we give the proof of Theorem [2l In 
Section [3] we apply inequality (fT6|) to derive the stability of the inverse problem presented in 
Theorem[T] The proof of Theorem[2]is organized in several subsections. In Proposition [T] we 
prove the properties of (j) which will allow us to use it as a weight in a Carleman estimate. 
In subsection 12.21 with 4> as the weight function, we develop the L^-product {Pi{w) , P2{w)) 
for functions with non-zero boundary values. We prove inequality (jl5|) in subsection (j2.3p . 

2 Proof of the Carleman inequality 
2.1 Weight function 

Here we prove that the function cf) satisfies enough properties for being a weight function in 
a Carleman estimate. 

We will use the following notation: 

M = Afilg^ +M2IQ2, a = oalQi +ailQ2, c{x) 



Qo^^ox (-T, T), = n,, X (-r, r). 

Proposition 1 // fii is a strictly convex domain of class , we can take e,6 > such 
that: 

(a) |V0| >6>0 m = {^^ U ^2) \ B,{xo) x (-T, T) 

(b) VMx,t)-Mx)>5>0, V(a;,t)eEi 

where vi is the unit outward normal vector to Qi. 
If additionally (I j jp is satisfied, we also have: 

(c) 4>i {x,t)^ 4>2 {x,t)^a2- (3t^ + Ml V (x, i) e El 

Qai+a2 gai+a2 

(d) "i g^aig^a2 ^ °^ i9a;"iga;"^ for all {x,t) e T.i and ai, a2 e NU {0} 
with ai + a2 < 3. 

(e) A(/)(x, t) > 2c{x) W{x, t) e . 

(f) D^{(l}){X,X)>6i\X\^ inQ^„ VX e for some 5i>0. 
Proof: We have 

V(/) — 2c{x){x — Xo) + |x — a:opVc(x). 
By definition, p{x) (thus also c{x)) is constant in the direction of a; — xq. Therefore 

(x — xq) ■ Vc(x) — 



and 



|V0|2 = Ac^{x)\x-xo\^ + \x-xo\'^\Vc{x)[' 
> Ac^{x)\x - xo\'^ 



2 

,2 



and (a) is proved. 
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Now, it is clear that (j)i{x, t) — a2 — Pt^ + Mi for each (a;, i) e Ei, so Fi x {t} is a level 
curve of 0i(-, t) for each t € [— T, T]. Since 0i(a;, i) < 02 — /3i^ + Mi < <j)i{y, t) for any a; € fii 
and y G ri2, we have Vi^i = |V0i|i^i on Ei, and thus (a) imphes (&). 

By definition p{x) = |a; — a;o| for all a; G Fi, hence (c) is simply deduced from (fTTj) . 

Without lost of generality, we can take = 0. Writing p in polar coordinates, Fi can be 
parameterized by 

7(61) = {p{e) cos B, p{e) sin 6'). (17) 

Qai+a2 

and then p is a function, li D = with ai + a2 < 3, we get 

0x1 0x2 

/ Ix — XqP \ 

aiD(/)i(x,t) = aifla-D f j = a2D(f>2{x,t) 

for all (a;, i) g Ei and (d) is proved. 

The expression for the Hessian matrix of second derivatives in polar coordinates is 

D^i^)^QgH{<P)Qj 

where Qg is the rotation matrix by angle 9, and 



( §1± 1 ( 9^'t> - 19± 

dr^ r \ drd9 r dO 



Now, we have that (recall that xo — 0) 



a 



One can notice that is well defined and smooth in ^Ixq, (which means {r > e} \ Fi). All 
the computations that follows are valid in this set. We already said above that p is constant 
with respect to r and only depends on 9 such that ^ = 0. Hence, we have that 



2a/ 1 -f 



H{<t>) = - ( M w,.2 , .2^ ) , (18) 



where we have denoted pg — ^ and so on. 

We will use the following well known facts (see for example |18j ) concerning curves in the 
plane: 

Lemma 1 Let 7 6e a curve in the plane. Then: 
(a) The curve 7 is strictly convex only at those points where their curvature is positive, 
(h) 7/7 is parameterized in polar coordinates by its angle, that is 

7(61) = {r{e) cos 9, r (9) sin 9), 
then the curvature of 7 is given by the formula 

r'^ -\- '2r'^ — rrgg 

^'>^^^^ = (^2+^2)3/2 • 
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Since the polar parametrization of Fi is precisely the above function with r{9) ~ p{6) and 
fii is strictly convex, we obtain 

^rAO)^'^0^>O We[0,2.[. (19) 

From HH) and HI]) we have 

A(f> = tr{D^{(t>)) 



tr(i/(0)) 

1 



and (e) is proved. 
We also have 



det(i?(</))) = ^{3pl~ppee + p 



2 



p2 

By the Sylvester's Criterion we can see that H{(f)) (thus D'^{(f))) is positive definite. Indeed, 
the element || and the determinant of the matrix i?(0) are positive. Finally, since fl is 
compact, this implies (/) and the proof of Lemma [1] is complete. ■ 

We introduce now the last hypothesis we will need in order to get the Carleman inequality: 

f minjai , flaWi Ml , ^ 

/3<min LAl^,_1 (20) 



and we take 7 g (0, 1) such that 



7 > . , '^ia. (21) 



/3- 



diam(n) 



2min{ai,a2}(5i 

^ 2/3 + max{ai,a2}||A</.|||^(^^/ ^ > 

Remark 2 1. We have to take M large enough in order to ^20^ and the hypothesis of the 
inverse problem (see Theorem\3j) become compatible. 

2. Taking [3 small enough, \21\ and l\22i become compatible. But, as we will see in the 
next section, the smaller is [3, the bigger need to be the inversion time for the inverse 
problem (see Theorem\B^. 

3. Actually, the optimal 5i is the first eigenvalue of (0) . Having an explicit expression 
for it could help to a better choice of (3. 

4-. From the hypothesis of Theorem\^we have 02 < ai and the maximum and the minimum 
in (1 1^01) and l\22i are known. 



Remark 3 Pseudoconvexity. It is easy to check that Proposition [7] and assumption 
imply that (f>i and and (j)2 are pseudoconvex fJl \2lV with respect to P in B^{xf)) and VL2 
respectively. Global Carleman estimates ( without explicit dependence on the parameter A ) can 
be deduced in each subdomain VLi \ Bg{xo) and il,2, for solutions that vanish on the exterior 
boundary and the interface (see Nevertheless, the traces of the transmission wave 

equation do not vanish on the interface Fi and we will need to use carefully hypothesis ai > 02 
and the parameter s. On the other hand, we have to use the parameter A in order to get rid 
of the lack of pseudoconvexity in B^{xo). 
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Remark 4 Coefficient a — a{x) variable. We can also prove a Carleman estimate in the 
more general case a{x) — ai{x)tni + a2{x)tn2 with Oj e C^{flj), j — 1,2, if each aj is 
constant in the interface Ti, and under some hypothesis on Va (similar to those of J2S^). 
More precisely, if we check the pseudoconvexity condition in this case, we will have that 4> is 
pseudoconvex with respect to the operator du — diviciS/u) (in each domain Qj) if there exists 
9 e (0, 1) such that: 



2aj 



<^i(l-e), j = l,2 (23) 

y/^j+T\Waj\ \ 2 J 

It is easy to check that these hypothesis are compatible with the assumptions of Theorem]^ 
Indeed, for T sufficiently large, there exists (3 satisfying both (3 > ai/T^ and {2^^ . Never- 
theless, in order to construct the weight function as we have done above ( and deal with the 
traces of the solutions on the interface Ti), it is crucial for each function a.j to be constant 
on the interface. 

2.2 Listing all the terms 

In this part of the work we develop the L^-product of Pi{w) and P2{vo). We will do formal 
computations, by writing generically (p for the weight function and Q for the domain with 
boundary S. 

As presented in Section 1, we have, for A > 0, s > 0, 

if^e^"^, w = e'^fu 



where 



P{vj) = e'*'^L(e^"'^w) ^ Pi{w) + Pziw) + i?(w) 
Pi{w) = wtt - aAw + s'^X^ip'^E{(j))w, 

P2{w) ~ {'J — l)sXLpL{(p)w — sX'^(fE{(j))w, —2sX(p{(j)tWt — aW(j) ■ Ww). 
R{w) — —^s\ipL{(j))w 

3 

We set {Pi{w) , P2{w)) = ^i,3^ where li^j is the integral of the product of the ith- 
term in Pi{w) and the jth-term in P2{w). Therefore, 

h,i = -sxh-i) ff ^miwti^ + '^'^Y^'^ ff \w\'^{^u + M<i>t\^)m, 



Iq J Jq 



bsX^ 
2 
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JJq JJq 

+sX // |w(p(/?a(A(;i + A|V0p) - sA // \wt\'^a(pV(f) ■ v, 
JJq J Jt, 

l2,i = -sA(7 - 1) JJ^ ipL{<j>)wa^ + sA(7 - 1) j a\VwWL{<l>) 

-sA2(7-1) jj \w\^aip{V^-VL{^)) + sX^^^ J J \w\'^aipVL{4)) ■ u 



'Q 

+sA2l^ // \w\^aipL{(j)) 



h,2 = sX^ 1 1 aipE{(j))wVw-v-^ // a\w\'^{X^pE{<i))V<i) + ^pWE{(i))) ■ v 



+^jj \w\''aipE{cp){A(P + X\V<P\^)-2sX^ JJ |w;|2aV^'(0)(V0, V</>) 
+^JJ \w\^a^A{E{cP))-sX^JJ \Vw\\^E{cl>), 



hfi = sxff \Vw\^a<pL{(j)) + sX^ [ [ \Vw\^aipE{(j)) +2sX^ I j aV|V(^-Vu;|^ 
JJq JJq JJq 

-2sA^ // aipcptWtVw ■'V<l) + 2sX // a'^<fiD'^{(j)){Vw,'Vw) 
JJq JJq 

+^^JJ \^wfa'^ipV(j)-i' + 2sX JJ aip{(j)twt- aV(j)-Vw)^, 
h,i = s'X^^-l)JJ^\wW'L{ct>)E{ct>), 
h,2 = -s^X" JJ^\w\\^E{<i>f, 

„3\3 // L,.|2, 3 El/ .\ r / -\ , o„3\3 // L..|2, 3/ij, |2 / , „2n2/ 



^3,3 = s'^A-^// \w\'^-'E{ct>)L{ct>)+2s^X' |u.| V^(l0t|>« + a •D"(</')(V<^, V(/.)) 
Iq JJq 



+3s=^A* JJ^ kl V'i^W' + s'A^ JJ^ a\wW^E{<l>)^. 
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Gathering all these terms, we get 

(PiH,P2H)i2(Q) = 2sxJJ\wt\^^(t,u-jsxJJ\ Wt\'^(pL{(l)) 

+2s\^ J I if {\wt\^\<j)t\^ - 2wt<j)ta\7w ■ \7(f> + a^\\7(f> ■ \7w\'^) 
+2sX J J a^(pD^{(t)){Vwyw) 



+2s3a3 // |«;|V'(|'/'t|'0tt+a'-D'(</')(V</>,V0)) 



+7s3A3yy \w\'ip^Licf>)E{^) 
+X + J, 

where J is the sum of all the boundary terms: 

J = sX ff (a^^\Ww\^^-2a^^{W<j>-Ww)^ 
JJt. \ ov dv 

+sA^^^ / / \w\'^ipa\/L{<P) ■ V 



2 

+sX'1—^ II \w\''ipaL{cP) 



s 

27-1 ff 2 



E 



-sA(7-l) // wa—(pL{<f)) + sX^ aipE{(t,)w— 



+2sA JJ atpcjjtwt^ - sX JJ \wt\'^(pa^ 



and X is a the sum of the remaining terms, in such a way that; 

\X\ < CsX"^ JJ ifi^lwf 

In the sequel, we denote by Aj, j — 1, ...,8 the first eight integrals we have listed in the 
product of Pi{'w) by P2{w). Thus, we have 

8 

(PiH,P2H>i2(Q) + J. (25) 

2.3 Proof of Theorem El 

We take /3, 7 and M satisfying the hypothesis pT|l . (1^ . (pij) and (1^ . and 4> as the corre- 
sponding weight function. We assume throughout all this part of the work the hypothesis of 
Proposition [T] (especially that Vli is strictly convex with boundary). 
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Recall the notation Qo = (^^i U l^a) x {-T,T) and = {VLq \ B,{xo)) x {-T,T). We 
apply the above computations to w = e^'^u with it G X in each one of the open sets Qi and 
Q2- Adding the terms that result in both cases, we have (recall that u{zLT) = ut{±T) — 
for all u e X ): 

8 

{Pi{w),P2{w))^,^Q^) ^Y.'^J'QoM + ^QoH + J^.H + Mw), (26) 

where we have written ^j,Qo instead of the integral Aj given in subsection 12.21 taken in the 
set Qo, etcetera. 

The proof of Theorem [2] is based on the next three facts: 

• The sum of the Aj-integrals in Q^^ can be minored. 

• The sum of terms in the interface given by Jsj^ is nonnegative, and: 

• We can introduce a second weight function centered at a point different to xq in order 
to deal with the integrals in B^{xq). 

The key points in each step of the proof are based on the properties of 4> listed in Propo- 
sition [TJ 

2.3.1 The interior 

Proposition 2 There exist 6 > 0, C > and Aq > such that: 

8 

' A - „ (n„\ ^ AILnll„ _ r'lUnll „ , , 



for all A > Ao, for all u Cz X . 

Proof: We arrange the terms into four groups: 

1. Ai^Qg + A2^Qo = sX |wt|V(-7-^(0) - 4/3). For all (x,t) in Q^^ we have: 



-7L(0)-4/3 = 7(2/3 + aA(?!)) - 4/3 

> 7(2/3 + 2ac{x)) - 4/3 by Proposition [1] 

> 7(2/3 + 2a-— -^^-—) - 4/3 by definition of c{x) 
= d>0 by 



Therefore: 



2. A3,Q„ = 2sA^ / / fi<j)twt " aV<j) ■ VwY > 
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3. ^4,Qxo+^5,Qx„ =sA // ip{a^D^{(j))(yw,Vw) + -faL{(j))\Vw\^) . Then, using Propo- 
sition [1] and (|^^ , we obtain 



> sA / / ^a(2a(Si -7(2/3 + aA(/))) |Vwp 



> sA// ¥3a(2a(5i -7(2/3 + all A0||l~)) iVwp 

> sAJ2 / / y^lVwp 



Therefore 



^4,Qo + ^5,Qo > '52SA / / ip\Vw\'-CsX ip\Vw\ 



4- E ^j.Qo = / / |«;|V'-Fa(0) where 

j=6 J JQo 



= 2XE{(j)f + -iL{(j))E{(j)) - IQP^t^ + 2a^D'^{(j)){V(l), 

= 2XE{(j)f + {iL{(j)) ~ 4:P) E{(t)) - 4/3a|V0p + 2a^ j:>^(0)(V(/), V0) 

bix)<0 

From Proposition [1] and (|20|) . there exists do > such that for all (a;, t) £ Q^g we have 

Fxic/)) > 2XE{(j)f +b{x)E{(j)) ~4:Pa\V(l)\^ + 2a^di\V(l)\^ 

> 2XE{(j))^ + b{x)E{(j)) + a{2a5i - 4/3)|V0p 

> 2XE{(^f - ||&||oo|i;(0)| + a{2a5i ~ 4/3)|V0|2 

> 2A^(0)2 - ||&||^|i;(0)| +do 
" h,x{E{<i>)) if£;(0) >0 

where 

X I — > 2Aa;^ + (-!)■' ||6||ooa; + do. 
As do > 0, there exists Ao > such that for all A > Ao 

mm(/,, > ^ > J = 1,2. 

Thus, for each A > Ao we have 



J=6 



E ^..Qo > ^3S^A3 / / kl V' - Cs'X' / / kl V' 



-8,(2:0) 



By collecting all the terms Aj^q^ together, we conclude the proof of Proposition [2l 
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2.3.2 The interface 

Since the interface Fi is a common boundary of fli and 172, the term J^^ in (|26p is the sum of 
the integrals comming from each domain: J^^ — J^iiwi) + Jy.i{w2)- We have the foUowing 
resuh: 

Proposition 3 Suppose < a2 < ai- Then there exists sq > such that 
Jsi = + Jsi(w2) > Vs > So 

/or all u X . 

Proof: We enumerate the ten integrals arising in (|25p associated with the common boundary 
Si, and we denote by Ji the sum of the i-th integral in ([25|l which comes from fii with the 
respective one of rt2- Ji — Ji{wi) + Ji{w2)- 

In order to prove the inequality, we arrange the terms into three groups. In each case, we 
use Proposition [T] and the fact that w satisfies the transmission conditions. 

1. Is not difficult to see that Jk = for each k £ {3,4,7,8,9}. Indeed, from (d) of 
Proposition [1] we get Li{(f>i) — ^2(02) and aiS/ Ei{(f>i) — a2Vi?2(02) on Ei, and the 
desired result follows. 

Now, let us denote by g the real function defined in Ei by 
gix,t) -.^ Ei{dpi) ~ E2{h) = 

Since 02 < ai, we have g > in Si. 
Thus we can prove: 

2. J2 + Je + 1^10 > Vs > So 
Indeed: 

—J2 ~ Je — s\ — ^— ^ J J \w\'^ ipaV L((j}) ■ v 




for all s > So, since Lp>l. 

3. J1 + J5 + 5J10 > Vs > So 

By construction, </> is constant on each level Fi x {i} of the interface (Proposition [T]) . 
Thus 

V0, .Vzz;, = ^5^ in Si for j = 1,2. (27) 
Moreover, since w satisfies ([5|) we have 



dwi 




dw2 






dT2 



in Si Vu G X. 



Hence: 

2 

E 




dTi 



dwi 



1 "i'^ ) ('^i - 02) > 



(28) 
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From ^ and ^ we get: 
Ji > -sA 
= -sA 
= sA 

On the other hand, 

-J5 = 



Si 



dwi 


' 901 ^ 


dw2 




<Pi^ 1- 


dv2 



V2 



dv2 



ai 



dwi 



ai 



dwi 



1 1 



(pi «i 

a2 ai ' 



901 
dvi 



1 



9wi A 
u^i^Uai— j 

1 





dwi 







Si 



ai 



dwi 
dvi 

\w\^ipig ( ai 



1 1 

(^1 ( 

02 ai 



Si 



90iA 



ai 



9wi 
dvi 



•^ly- ^ ) 1 ai 

a2 ai 



1 



< -j^Jw + Ji Vs > si. 



Proposition Q is proved. 



2.3.3 The boundary S. 

Since we deal with functions w such that W2 = in E, we have 



Jt. = Ji[w2) 



0-2 



dw 



> -sA 

> -sA 



dw 




a2^ 
av 


\ 


dw 


2 


dv 





dv 

d. 
dv 

d(t> 



dv 



L~(S) 



-sAC 



0.2 



dw 
dv 



where we have defined S+ = {(x, t) G F : V(/)(x, t) ■ v{x) > 0}. 
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2.3.4 Carrying all together. 

From (|26p . pS)) and Propositions [5] and [31 there exist sq, Aq, C G M such that for each s > sq 
and A > Aq we have 



M -c\\wf 



sAC 



02 



9:/ 



< 



Adding ^ I Pi HI 



2 



at both sides of ([50)1 we obtain 



i^iHIl^(qo) + i^2Mi,.(q„) + iiz«ii;,^-c 





dw 




Of 





<c|PiH + F2H|^.(Q„) 
<c(|PHli.(Q„) + |PHli.(Q,)) 

Thanks to ([20]) we have A < C(p for A large enough. Therefore 

<^||w||^ Vs>si. 
From (|3ip and ([5^ we get, for all s > max{so, si}, A > max{Ao, Ai}: 



\PlH\hiQo) + \P2H\lHQo) 



W 



+^IIHlL.„,,,+-AC//^^^|a2|^| 



2.3.5 Eliminating the term in B^{xo). 



(30) 



(31) 



(32) 



(33) 



In the last step we will remove the integral in Pe(xo) from the right hand side of (|33p . In 
order to do that, first remark that xq can be arbitrarily chosen in lli since fii is strictly 
convex. 

Thus, we can take two different points in Qi and we have the two respective inequalities 
given by ((33)) . Now, we will show that the left hand side of each inequality can absorb the 
term || • ||_Be(a;o) from the other inequality provided that e is small and A is large enough: 

Denote by xi, X2 two points in f2i, and <p^, (jp' their respective weight functions. In order 



to have || • absorbed by the term 



it suffices that 



in B^{xi) 



gAW=-0i) ^ 2C in B,{xi) 

Thus, if we show that it is possible to have (fp' — (f)^ > 5 > Q \n B^{xi) by taking A large 
enough we are done. 

In fact, let he d — ^\xi — X2\ and assume that e < d. Then, for all a; G Bs;{xi) we have: 



<j)\x,t) < —e^-(3f + M 
Pi 

< - pt^ + M, 



(34) 
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where ai = d{xi,Ti) > 0. 

In the same way, if we denote D2 — max d{y, X2), we get 



> 



Consequently, we have 



(3f 



M 



Pt^ +M 



> a 



^2 



Vx e B^{xi). 



(35) 



It is clear that an analogous result is true by interchanging xi and X2 (now with a2 and 

dai da': 
D2 ' Di 



Theorem [2] is proved. 



Di). Thus, taking e < min ( ^yJ-, j -^^g ^^n absorb the desired terms in the inequality and 



3 Proof of the stability of the inverse problem 

In this section we apply the Carleman inequality of Theorem [1] to the inverse problem pre- 
sented in Section [TJ For a principal coefficient a piecewise constant and p G L°°{n), we 
consider the wave equation 



Utt — div(a(x)Vu) + p{x)u 

u 

uiO) 
ut{0) 



g{x,t) 
h 

Ul 



n X (o,T) 

Fx (0,T) 

n. 



(36) 



If .g e L\0,T;L^{n)), h £ L^{0,T; L^{T)) and uq £ H^i^), wi £ L^{^), then [32l[l3| equa- 
tion ((36|) has a unique weak solution u e C{[0,T]; H^{n))r]C'^{[0,T]; L'^{n)) with continuous 
dependence in initial conditions and such that ^ £ L^(0, T; L^(F)). 

In order to prove the local stability of the nonlinear application ([2|), that is, the problem 
of determining the potential p in by a single measurement of the flux a2^ on F between 
t = and t ~ T, wc follow the ideas of 9J and [27j Thus, we will first consider a linearized 
version of this problem, what means working on the wave equation 



ytt - div{a{x)Vy) 



-p{x)y = fix)R{x,t) nx{0,T) 

y = rx(0,T) 

2/(0) = n 

yt{0) = n 



(37) 



given p and R, and proving the stability of the application /|o 
We will indeed prove the following result: 



Theorem 3 For xi,X2 £ fli let Rj = sup{|a; — yj{x)\ : x £ ft2}, j — 1,2, whe 



with xq = Xj . Set aj 



d{xj,Ti) and Dq — max 



defined in 

- — ■ \ cti a2 

the hypothesis 0/ Theorem]^ and T , (3 satisfying and suppose that 



Ri + ai R2 + a2 



With 



• IIpIIl-(o) < m 

• R£ H\0,T;L°°{n)) 
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• < r < 0)1 almost everywhere in fl. 

Then there exists C > such that for all f G L'^(fl), the solution y of ( | j*?] ) satisfie 

2 



li2(o) < 



5y 



Proof: For each / e L2(f]) and i? £ i?i(0, T; L~(fJ)), let ?/ be the solution of We take 
the even extension of R and y to the interval (— T, T). We call this functions in the same way, 
and in this proof we denote Q — [~T,T) and E = F x {—T,T) the extended domains. 
Therefore, z = yt satisfies the following equation: 



ztt - div(aVz) + = f(x)Rt{x,t) Q 
z = E 

z(o) = n 

zt{0) = f{x)R{x,0) n 
and we have the usual energy estimate 

\\z\\H^-T,T-L^{n)) < C\\f Rt\\L^{^T.T;L^{ny) + C||/-R(0)IIl2(O) 

that gives, since R e H\0,T; L°^{n j), 

\\z\\H^(-T,T-L^in)) < C\\f\\L2(n)\\R\\H^{0,T:L°°{a)) < ^11/11^2(0). 



(38) 



(39) 



In order to apply Theorem [2] and use the appropriate Carleman estimate, we need a 
solution of the wave equation that vanishes at time t — ±T. Thus, for Q < S < T we take 
the cut-off function e Cg°{-T,T) such that 



• < 6* < 1 

• 9{t) = 1, for all t e (-T + S,T - 5) 
and we define v — 9z. Then v satisfies: 



' vtt - div(aVt.) + = e{t)f{x)Rtix,t) + 2etytt + 0ttyt Q 

V ^ S 

v{o) = n 

vtiO) = f{x)Rix,0) n 

v{±T) = vt{±T) = n. 



(40) 



Take j e {1,2}, and let y be the function defined in ^ and (j) the weight function, corre- 
sponding to the point Xj G fli. Notice that 

0(x, t) < (j){x, 0) V(a;, t) e (0, T) x fl. (41) 

Moreover, by definition of p and a (see ^ and the definitions below) we also have 



p{x) 



< 1 



\x - yix)\ 
p{x) 



< 1 



< Do 



and then 



(l){x,t) 



(3t'^ + M 



< aDl - I3t^ + M. 
Then, by the choice of T > Dq^J^ we get 

(j){x, ±T) <M < 0(x, 0). yxen 



(42) 



(43) 
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Thus, taking 6 small enough, it is also true that 

^{x,t) < M <(l3{x,0). (44) 

for allxefl and t e [-T, -T + S]U[T- 6, T]. 

From now on, C > will denote a generic constant depending on f2, T, 9, xi, X2, S, 
So and Aq but independent of s > sq and A > Aq. We will occasionally use the notation dt 
for the time derivative. 

As in the proof of Theorem [21 we set ip — e^'^, Wj = e'"''wj and 

Pi{w) = wtt - aAw + s^X^Lp^E{(j))w. 

It is easy to check that 

2 



_ 1 



\dtw{fd)\'' + X 



(45) 



where X is a sum of negligible terms such that 

1 



Jn 



Since we have u;t(0) = e'*'^(°)wt(0) = e'"^(o)/(x)i?(a;, 0) and |i?(x,0)| > r, we get (recall that 
Q = f7 X (-r,r), = % X (-r,r) and so on). 



In order to apply the Carleman estimate (Corollary [2]) we consider both weight functions ip^ 
and <p^, corresponding to Xi and X2 € ^^i and we apply the previous estimates to = e^^ Vj 
for j, k = 1,2 and sum up the inequalities. We obtain, for s > sq and A > Aq, using Cauchy- 
Schwarz inequality, the following: 

^2 f (g2V(0)+g2V(0))|J|2 



< 



fc=i v=i 



ytw. 



L2(njX(0,T)) 



fc=l 



< 



^ i,fe=l 

Now, applying Corollary [21 we get 



i.||2 



r' (e 



c 



C 



fc=i 

2 



^ fc=i ^ 



fc^3|„,,fc|2 



(^'^)1^ 



02- 
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On the one hand, since we have 6*4 = in [— T + S,T — S], then from estimate (l39l) . (|4T|) and 
(I44l). we obtain Vfc= 1,2 



Qo 



< c 



Qo 



Qo 



Qo 



„2V(0)|J|2|^^|2^C' 



-T+S 



< 



T-s J Jno 

Jn 

Jn 



{\^t\' + \z\' 



T;L2(n)) 



Now, recalhng the notation stated in p^ . we have U S!^ C S and 
S for each k, so we finaUy obtain 



dw'' 




dv 


\ dv \ 



^2 / (g25yi(0) g2sv^(0))|J-|2 

C 



< 



(0) _^g2V(0) 



I/I' 



(46) 



dv 



For s large enough, the left hand side in can absorb the first term of the right hand 
side. Therefore, since Lp^ and 9 are bounded on S and ^ is an even function with respect 
to t g [-r,r], we obtain 



I/I' < C 



< C 



= 2C 



/ 1 2s(fl^ I 2 2s(/3' 



dv 
av 



a2- 



dz 







and this ends the proof of Theorem [3l ■ 

We will end this paper by the proof of Theorem [1] which is a direct consequence of 
TheoremO Indeed, if we set y — u{q) — u{p), f = p — q and R ~ u(p), then y is the solution 
of 

yu-div{a\/y) + {p- f)y = f{x)R{x,t) {Q,T) x n 

y = (0,T)xE 

y(0) = O ^ '> 

ytiO) = 17 

where q — p — f E U, with bounded in L°°(ri) from the hypothesis of Theorem[T] The key 
point is that in the proof of Theorem [3l all the constants C > depend on the L°°-norm 
of the potential as stated in Corollary O Thus, with q G U, we are actually, with equation 
(|47p , in a situation similar to the linear inverse problem related to equation (j37p and we then 
obtain the desired result. ■ 
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